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Partial widths of the radiative decays 0(1020) 7/0 (980), 777, 777', 'j-k'^ and /o(980) — + 77 
are calculated assuming all mesons under consideration to be qq states: (/)(1020) is dominantly 
an ss state {nn component < 1%), J), rj' and vr" are standard qq states, rj = nn cos 9 — ss sin ^ 
and rj' = nn sin 6 + ss cos with 9 ~ 37°, and /o(980) is the qq meson with the flavour wave 
function nncosip + sssin^. Calculated partial widths for the decays 0(1020) — > 77;, 777', 77r'' are 
in a reasonable agreement with experiment. The measured value of the branching ratio BR{cf> — > 
7/o(980)) requires 25° < \(fi\ < 90°; for the decay /o(980) — > 77 the agreement with data is 
f~>) ■ reached at either 77° < ifi < 93° or ( — 54°) < (p < (—38°). Simultaneous analysis of the decays 

' 0(1020) ^ 7/o(980) and /o(980) — » 77 provides arguments in favour of the solution with negative 

. mixing angle ip = —48° ± 6°. 
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The problem of the /o(980) structure is actively discussed during the last decade, see for example and references 
therein. The question is whether /o(980) belongs to the scalar quark nonet l^Poqq or whether it should be considered 
, as an exotic state. Recent measurments of the radiative decay 0(1020) 7/0 (980) [0,D have reinforced the discussion. 
As was stressed in pP,p^, it is reasonable to perform the qq nonet classification in terms of the matrix, or 
"bare", poles corresponding to states "before" the mixing which is caused by the transitions qq state — > real mesons. 
^""^ Such a mixing is crucial for the formation of the scalar/isoscalar resonances in the region 1200-1600 MeV, which 
^ are /o(1300), /o(1500) and the broad state /o(1530l;25o)- According to the X-matrix analysis, these resonances are 
formed as a result of mixing of the lightest scalar gluonium with the l^Poqq and 2^Poqq states. 
^ ; In terms of f^"''^, the scalar/isoscalar states of the basic l^Poqq nonet are /o^"'''=(720± 100) and /o^'''''=(1260±30); the 
^pH, state /q'"''^ (720) is close to the flavour octet ip — —69° ±12° (the mixing angle is determined as /o = nn cos ip + ss sin ip 
O 1' where nn — {uu + dd)/\/2), while the state /g'"''^(1260) is almost a singlet (ip — 21° ± 12°). The transitions which are 
D - due to the decay processes /o tttt, KK, r/rj mix the states /q°'''^(720) and /q'"''^(1260) with each other as well as with 
nearby states: jfo^'^''^(1235±50), /o^''''^(1600± 50), /o^''''^(1810± 30) (the orthoganality of the coordinate wave functions 
^ I does not work in the transition {qq — state)i — > real mesons — > {qq— state)2) ■ Two of these states are the 2^Po'?9 nonet 
members and one is the scalar gluonimn, either /q"''^(1235) or /q°'''^(1600): the iiT-matrix analysis cannot definitely 
. , tell us which one. Lattice calculations ijllj favour /q'"''^(1600) as a candidate for the glueball. After the mixing, 
' we have in the wave {IJ^^ = 00++) four comparatively narrow resonances (/o(980), /o(1300), /o(1500), /o(1750)) 
" ■ ' and one rather broad state /o(1530^25o)- '^'^^ broad state (r/2 ~ 400 — 500 MeV) appears as a result of a specific 
phenomenon which is appropriate for scalar/isoscalar states in the region 1200-1600 MeV, that is the accumulation of 
widths of overlapping resonances by one of them ]l^ . The analogous phenomenon in nuclear physics was discussed in 
p3| . Fitting to data proves that the main participants of the mixing are the glueball and two qq states belonging to the 
nonets l^Poqq and 2^PQqq, whose flavour wave functions are close to singlets: the transition glueball —> qq{singlet) 
is allowed, while the transition glueball — > qq{octet) is suppressed. Gluonium/gg mixing is not suppressed according 
to the 1/A^ expansion rules JT^ (see |^,^| for details). The glueball descendant is a broad state /o(1530l^^o) which 
contains about 40-50% of the gluonium component; the rest is shared mainly by /o(1300) and /o(1500). 

Following the ii'-matrix scenario, the resonance /o(980) is the descendant of /p*'"'^ (720 ± 100): the shift of the bare 
pole to the region of 1000 MeV is due to transitions of /q'"''^(720 ± 100) into real two-meson states, tttt and KK. 
As to its origin, the state /o(980) is the superposition of states qq, qqqq, KK and tttt, and the wave function is to 
be represented by the correspondent Fock column. Concerning the KK and tttt components, one should take into 
account only the part of wave functions which respond to large meson- meson separations, r meson > Rconfinement] 
on the contrary, for the qqqq state small interquark separations are to be taken into account. As to the glueball 
admixture in /o(980), the relative suppression of the decay J /'^ ^ 7/0 (980) tells us that it is not large. 

To enlighten the problem of the content of /o(980), it would be reasonable to suggest a dominance of certain 
component: basing on the results of the i^T- matrix analysis P,pyiC|] it is natural to check whether the hypothesis about 
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the dominance of qq component can explain the data on radiative decays (/)(1020) — > 7/o(980) and /o(980) 77. 
This study is dictated by the existing opinion that measured widths of radiative decays contradict the interpretation 
of /o(980) as qq state, see mini-review and references therein. 

In Section 2 we calculate the decay amplitude 0(1020) 7/0 (980) assuming /o(980) to be a qq system. The 
technique for the description of composite qq systems and corresponding form factor calculations has been presented 
in detail in Ref. ||l^, where the pion form factor was studied as well as transition form factors tt" 77*(Q^)i 
ij — > 77* (Q^) and 77' 77* (Q^)- The method of spectral integration has been used for the form factor calculations. 
Within this technique, it is possible to introduce normalized wave functions. The transition amplitudes obey the 
requirement of analyticity, causality and gauge invariance. The used techique allows us also to perform calculations 
in terms of the light-cone variables. 

It is worth noting that this calculation technique for the processes involving bound states has a broader applicability 
than for qq systems only: in [ p^ this method got its approbation by describing the deuteron as a composite np 
system, then this very technique was applied to heavy mesons p0| , |2l[ |. For the convenience of a reader, in Section 
2 we recall briefly the basic points of this approach and give necessary formulae for the calculation of the reaction 
^1020) 7/o(980). Our calculations show that the data on branching ratio BR{<j) 7/o(980)) = (3.4±0.4to 5)-10-'' 
[^,p|Jl^ may be described within the qq structure of /o(980), with the mixing of ss and nn components as follows: 

ipfiavour{fo{Q80)) = nn cos (fi + ss sin ip , 25° < |<^| < 90° . (1) 

The width r(0 7/o(980)) depends strongly on the radius squared of /o(980), ^^^(ggo)- At large -^^^^(ggo) > 12 GeV^^ 
(0.47 fm^), the data require \ip\ - 25° - 45°, while at i?y„(98o) ~ 8 GeV-^ (0.32 fm^) one needs \(p\ ~ 40° - 75°. 

The estimation of radii of the scalar/isoscalar mesons, which has been carried out by using GAMS data for tt^p 
7r°7r°n shows that qq component in /o(980) is comparatively compact, -R^^jggo) = 6 ± 6 GeV~^ Q. So, basing 
on GAMS data, it is reasonable to put 

6 GeV-2 < Rl^^^^^ < 12 GeV-2 , (2) 

that favours the large values for the mixing angle. 

Concerning the scheme under investigation, it is of principal importance to verify if the transition (/)(1020) — > 
7/o(980) — > 77r°7r'' describes the measured tt^tt^ spectrum within the Flatte formula with parameters determined 
by hadronic reactions. In Section 3 we calculate the tt^tt" spectum using the parametrization of the Flatte formula 
suggested in |^,^, and reasonable agreement with data is obtained. The performed spectrum calculations indicate 
the existence of significant systematic errors in extracting the resonance signal 0(1020) — > 7/0 (980) which are caused 
by the background contribution. 

Radiative decay /o(980) — > 77 provides us another source of information about the content of /o(980): the cal- 
culation of partial width /o(980) — > 77 has been performed in |2^ assuming the qq structure of /o(980). However, 
as was stressed in Q , the extraction of the signal /o (980) 77 from the measured spectra 77 tttt meets with 
a difficulty, that is, a strong interference of the resonance with background, thus giving uncontrollable errors. The 
recently obtained partial width r(/o(980) -f-f) = 0.28;[! ?^ keV [PStl is a factor two less than the averaged value 
reported previously (0.56 ± 0.11 keV [M)- In Section 4 we re-analyse the decay /o(980) — > 77 using the Boglione- 
Pennington number for the width (0.28^q'53 keV) together with the restriction for the mean radius squared of /o(980) 
p3| : ^^^(ggo) < 12 GeV^^. This provides two possible intervals for the nn/ ss mixing angle: 

/o(980) 77 : 1) 80° < < 93° , 2) (-54°) < Lp < (-42°) . (3) 

In Section 5 we calculate partial widths for decays 0(1020) — > 777, 777', 77r'^ in the framework of the same technique 
as is used for 0(1020) 7/0 (980) and with the same parametrization of the 0-meson wave function. The calculations 
demonstrate a good agreement with the data for these processes as well. The calculation of the decays 0(1020) — > 
777, 77;' provides us with a strong argument that the applied method for calculating radiative decays of qq mesons is 
wholly reliable, and the decay 0(1020) 77r° allows us to estimate the admixture of the nn component in 0(1020). 
The partial width of 0(1020) 7ao(980) is also proportional to the probability of the nn component in 0(1020), we 
discuss this decay in Section 5 as well. 

The restrictions for the mixing angle Lp which come from the combined analysis of radiative decays 0(1020) — > 
7/0 (980) and /o(980) — > 77 are discussed in Section 6. For negative mixing angles the combined analysis does not 
change the restriction (3): 

p = -48° ± 6° , (4) 
and does not provide any restriction for the radius -R^q(98o)('^"-)- 
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For positive mixing angles the allowed region of ((^, -R^Qjggo) ('^"')) reads as follows: 

^ = 86° ± 3° , i?^^^(ggo)(nn) < 6.8 GeV"^ . (5) 

This means that in this solution /o(980) should be rather compact ss state. 

The hadronic decays of /o(980) strongly contradict the almost complete absence of the nn component in /o(980), 
that gives the priority to the solution (4), according to which /o(980) is close to the octet state (recall, ipoctet = —54.7°). 

II. (^(1020) ^ 7/o(980): THE DECAY AMPLITUDE AND PARTIAL WIDTH 

In this Section we calculate the couphng constant for the decay (j) 7/0 (980) assuming the qq structure of /o(980). 
In the calculation we use the spectral integration technique developed in iQ . Within the framework of this approach, 
we study the reaction with a virtual photon, (p 7(g^)/o(980), and then the limit q^ ^ is attained. 

A. Transition amplitude 

The transition amplitude 4> ~^ 7('Z^)/o(980) contains S and D waves. Correspondingly, the spin-dependent part 
of the amplitude, (the index fi stands for the spin of the photon and v to that of the 0-meson), consists of two 
terms which are proportional to g^^^ {S wave) and q^q:^ — q±g'tu/-^ wave): 



(6) 



Here q'^ is the relative momentum of the produced particles (it is orthogonal to the meson momentum p) and gj^^ 
is a metric tensor in the space orthogonal to p: 

9^ = q^^ -P^^^ ' = - ■ (^) 

The requirement q^A^i, — results in the following constraint for the invariant amplitudes Fsiff') and Fu(q^): 

Fs{q') + l{qq^)FD{q') = , (8) 



hence equation (|^) can be re-written as 

A,. - leFsiq') (^g^. - ^ ^^0^7/0(9') 5^.^ • (9) 

The metric tensor g-^^ works in the space orthogonal to the reaction momenta: gj^^Pi, — and g^^q^ = 0. 

B. Partial width 

The decay partial width is determined as 

m^r^^/n^ = d^{p;q,kf)\A^,\^ . (10) 

In (0) the averaging over spin projections of the (p-meson and summing over photon ones are carried out; is the 
two-particle phase space which is defined as follows: 

For the radiative decay 4> 7 + /o(980), one has / d$(p; q, fc/) = (to^ — Mj) / {Idnm'^) . 
After summing over spin variables, the partial width reads: 

1 ml- Af? 

m^r^_f,^ = -a ^ , ^ \A^^^f,{())?. (12) 
D 

Here a = e^/An = 1/137. 
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C. Double spectral representation for the transition form factor (^(1020) 



7/0 (980) 



Assuming the qq structure of /o(980), the couphng constant of the decay 4> 7('7^)/o(980) is determined by the 
processes cf) ^ qq and qq /o(980), with the emission of a photon by the quark in the intermediate state, see Fig. 
la. Following the prescriptions of Ref. [Q, we calculate the three-point quark loop shown in Fig. lb by using the 
double spectral representation over intermediate qq states marked by dashed lines in Fig. lb. 

To be illustrative, let us start with the three-point Feynman diagram. 



1. Double spectral integral for three-point quark diagram 



For the process V jS, where V and S stand for scalar and vector particle, the Feynman diagram of Fig. la 
reads: 

AiFeynman) _ Jl_2_ r< ^V-^fS '^A"^ . . 

'^'^ J i(27r)4 [m? - kl){m^ - k'^){'m? - k^) ' ^ ' 

Here fci, k'l^ are quark momenta, m is the quark mass, and Gy, Gs are quark-meson vertices; the quark charges 
7(99) 



are included into Z^^^2-ys- The spin-dependent factor reads: 



4i:~^''^ = -Sp[(fcl+m)7^(fci+m)7,^(-A;2+m)] , (14) 

where the Dirac matrices 7^ and 7jj" are orthogonal to the emmitted momenta: ^jj^q^i — and ^^Pu — 0. 

To transform the Feynman integral ( p^ into the double spectral intcrgral over invariant qq masses squared, one 
should make the following steps: 

(i) consider the corresponding energy-off-shell diagram. Fig. lb, with = (fci + fc2)^ > 4m^, P'^ ~ {k[ + fc2)^ > 4m^ 
and fixed momentum transfer squared q^ — {P — P')^ , 

(ii) extract the invariant amplitude separating spin operators, 

(iii) calculate the discontinuities of the invariant amplitude over intermediate qq states marked in Fig. lb by dashed 
lines. 

The double discontinuity is the integrand of the spectral integral over P^ and P'^. Furthermore, we put the following 
notations: 

P2 = s, P'2 = s' . (15) 

For the calculation of discontinuity, by cutting Feynman diagram, the pole terms of the propagators are replaced with 
their residues: (m^ — fc^)~^ 5{m? — fc^). So, the particles in the intermediate states marked by dashed lines I and 
II in Fig. lb are mass-on-shell, fc^ = fc| = k'^ = m^. As a result, the Feynman diagram integration turns into the 
integration over phase spaces of the cut diagram states. Corresponding phase space for the three-point diagram reads: 

d$(P, P'- ki,k2,k[) = d$(P; kuk2)d<P{P'; k[,k'^){2nf2k2o5^'^\k'^ - k^) , (16) 

where the invariant two-particle phase space c?<I>(P; ki, /C2) is determined as follows: 

The last step is to single out the invariant component from the spin factor (O) . According to (|^) , the spin factor ( p^ ) 



is proportional to the metric tensor: Sj^^^^ ~ d^i!'- Then the spin factor Sy'_l^g, determined as 



is equal to: 



Sll-'^=9U^S^^X,{s,s',q% (18) 



S^^%{s, s', q') = -2m ( 4m^ + . - s' --^f_j«(s, s' , q' 
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2g^(s + s') — (s — s')^ — 

Here we have used that (fcifc2) = s/2 — m^, {k[k2) = s'/2 — m? and (fcj^fci) = — q^/2. 

The spectral integration is carried out over the energy squared of quarks in the intermediate states, s — — 
(fci + fc2)^ and s' = P'^ = {k[ + fc2)^, at fixed q^ = {P' — P)^. The spectral representation for the amplitude 

^{F„an) ^ ^j_^^^iFeynn^an)^^2^ rCads: 



OC 



ds f ds' Gv Gs 



Ay^^s - J ^ J ^ s- Mis'- Ml 

4m 2 4m ^ 



X j d<^{P, P'; fci, k[, k2) 4%s(s. 'Z')4'l^s • (20) 

It is reasonable to name the ratios Gy /{s ~ My) and Gs / (s' — Af|) the wave functions of vector and scalar particle, 
respectively: 



^=^v{s), -i%F2=^s{s'). (21) 



s-M^ ^ " s' -Mg 



Equation ( pO[ ) is a spectral representation of the Feynman amplitude shown in Fig. la for point-like vertices. In 
general case it is necessary to take account of the s-dependent vertices of the reaction. 



2. Form factor 0(1020) 7/0 (980) 

Generally, the energy-dependent vertices can be incorporated into spectral integrals. According to ||lj,|l^, the 
form factor of a composite system can be obtained by considering the two-particle partial-wave scattering amplitude 
1 -t- 2 1 + 2, with the same quantum numbers as for the composite system. In this amplitude, the composite system 
reveals itself as a pole. The amplitude for the emission of a photon by the two-particle-interaction system has two 
poles related to the states before and after electromagnetic interaction, and the two-pole residue of this amplitude 
provides us form factor of the composite system. When a partial-wave scattering amplitude is treated using the 
dispersion relation N / D-method, the vertex G{s) is determined by the A^-function: the vertex as well as A^-function 
have left-hand-side singularities which are due to forces between particles 1 + 2. 

Generally, the form factor for the tansition 0(1020) 7((7^)/o(980) reads: 

00 00 

^0-7/0(9^) J J -^V'0(s)V'/o(s') 

4m^ 4m^ 

X I d^P, P'- fci, k[M) Si%,is, s', q')Z^;%„ . (22) 

The factor z'^^Sfa defined by the fiavor meson component: for qq — nn we have -^^^/o ~ V2/3 and for the ss 

component ■^0'^/,, = —2/3. 

Working with Eqs. (p^), one can either 

(i) express it in terms of the light cone variables, or 

(ii) keep the spectral integrals over s and s' and eliminate integrations over quark momenta with the help of the phase 
space 5- functions. 

a) Light-cone variables 

The formulae (|l6|), (p7|), ( p2| ) allow one to make the transformation to the light-cone variables, using the boost 
along the z-axis. Let us use the frame where the initial vector meson moves along the z-axis with the momentum 
p 00: 

P={p+±,0,p), P'^(p+fl±^,_5-^,p). (23) 
In this frame the two-particle phase space is equal to 
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X 5 



and the phase space for the triangle diagram reads: 

1 dxidx2 



d<i>iP,P';ki,k2,k[) 



I6n x\x2 
X (5 I s — 



d^ki^d^k2±5{l ~xi- x2)5^^Hkii^ + k2A_) 



Xi 



X2 



m? + — (f±)^ + k? 



After having integrated over s and s', by using (5- functions, we have 

dx 



^0-7(9")/o('?' 



2^ _ ^0^7/o 



167r3 J x{l-xY 



X2 



(24) 



(25) 



(26) 



where x = k2zlp , k±_ = k2±, and the qq invariant energies squared are 



+ kf , + (k^ + xq±_) 

s = — -, =r, s = 



x{l — x) 



x{l — x) 



(27) 



b) Spectral integral representation 

In Eq. (|2^) one can integrate over the phase space keeping fixed the energies squared, s and s' . After using the 
phase space (5-functions, we have: 



dsds 



■V'0(s)V'/o(s') 



4m^ 



16\/A(s,s',Q2) 



with 



A(S, S', Q2) ^ (s' _ s)2 + 2Q2(5' + 5) + , 



(28) 



(29) 



The 6'-function, 0{X), restricts the integration region for different = —q^: d{X) = f at X > and 0{X) = at 
X < 0. 

In the hmit — s- 0, the integration over s' is carried out, and we have: 



ds 

^0^7/0(0)= / — ^0(s)V'/o(s) 



4m2 



■In- 



27r ^ - - 4m2 47r 



+ - — \/ s{s — Am?) 
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(30) 



To calculate ^^^^/^(O) one should determine the wave functions of the 0-meson and /o(980). 



3. Wave functions o/<^(1020) and /o(980). 

We write down the wave functions of 0(1020) and /o(980) expanding them in a series as follows: 

^ij,{s) = [nnsiwipv + ss cos 'pv)'^<t>{s) , 

^'/(s) — {nn cos (fi + sssiiup) ipf{s) . (31) 
For ipffi^s) and ipf{s) the exponential parametrization is used: 
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(32) 



The parameters and bf characterise the size of the system, they are related to the mean radii squared of these 
mesons, and ^jj^fggg) ^-nd and Cf are the normahzation constants. In our approach the dynamics of interaction 
is hidden in the wave functions and ipf^. We assume exponential s-dependence for the wave functions. Recall 
that exponential parametrization is often used in the quark model. Working upon meson radiative decays, we have 
additional arguments to employ this presentation of wave functions of the low-lying qq states. 

In Ref. [Q, in our analysis of pion form factor as well as transition form factors 7r° 11* {Q"^), V ~^ 11* iQ"^) 
and ij' 77* (Q^), we have reconstructed wave functions of pseudoscalar mesons as functions of s. At small relative 
momenta of quark and antiquark, < 0.6 (GeV/c)^ (where = s/4 — m^), the wave functions demonstrate rapid 
- exponential-type - decrease, then it becomes slower. The same behaviour is appropriate to wave function of the 
photon, which was used in our calculations; it is shown below, in Fig. 5). We keep in mind just this form of the 
wave functions for low- lying qq mesons, that is a foundation for the formula (^2|). In at calculation of the 
transitions r/ 77, ry' — > 77, ry(1295) 77 and ?7(1440) 77 we checked the possibility to use exponential-type 
wave functions. The processes r/ 77 and ry' — s- 77 have been calculated both with the wave function reconstructed 
from the experiment and exponential-type wave functions, with the same R^. The difference between form factors is 
not greater than 5% . 

At fixed and i?^^ ^ggp-j the constants and C / are determined by the wave function normalization, which itself 
is given by the meson form factor in the external field, Fmesoniq^), namely, at small q^: 



FmesoniQ ) — !-)- ^RmesonQ ■ 



(33) 



The requirement FmesoniO) — 1 fixes the constant 
parameter bmeson- 

The form factor Fmesoniq"^) reads |l8||: 



meson 



in (p^, while the value i?^ 



is directly related to the 



F„ 



dsds' 



e{ss'Q--m-X(s,s',Q-)) 



(34) 



where Smelon for /o(980) and 0(1020) is determined by the following traces: 

2P^sf;^ {s, s', q') = -Sp [Ck[ + m)7^(fci + m)(-fc2 + m) 
2P^^4*'-'(s,s',g2) ^ -isp [7;^^(fc; +m)7^(fci +m)7^(-^2 +m) 



(35) 



Here the orthogonal components are as follows: 



p 



1^ =7m-'?a'^ 

la — la ~ Pa j 7a la ~ Pa 



(36) 



and q = k[ — ki. When determining S^^^\ we have averaged over polarizations of the (j) meson (the factor 1/3). 
The functions S'j*'^'' and S'^^^^ are equal to: 



4*'-) (s, s', q') = ais, s', q'){s + s'- Sm^ - q') + q' , 
Si'"-^ {s, s', q') = I [a{s, s', q'){s + s' + 4m' q') + q'] 



(37) 



The radius squared of the nn component in the (/)-meson is suggested to be approximately the same as that of the 
pion: R^{nn) ~ 10.9 GeV~^, while the radius squared R^{ss) is slightly less, R'^{ss) ~ 9.3 GeV^^ (that corresponds 

to b^ = 2.5 GeV~^). As to /o(980), we vary its radius in the interval 6 GeV"^ < R%(^gso) < 18 GeV"^. 
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4. The result for (?i(1020) -> 7 + /o(980). 
The amplitude A^^^f^ (0) is a sum of two terms related to the nn and ss components: 

^0-7/0 (0) = cos (p sin (fivA'^^Zl (0) + sin (p cos <Pv^0^^/o (0) . (38) 

In our estimations we put cos(/?y 0.99 and, correspondingly, | sin(y5\/| < 0.1; for /o(980) we vary the mixing angle 
in the interval 0° < \ip\ < 90°. 

The results of the calculation are shown in Figs. 2 and 3. In Fig. 2 the values (0) and ^l^ll.yfg (0) are plotted 

with respect to the radius squared i?^^ ^gg^^ , while the mean radius squared of the 0- meson is close to that of the pion: 
^0(1020)- 0-4 fm2. 

In Fig. 3 one can see the value BR{(j) 7/o(980)) at various sin|(^|: sin|(p|=0.4, 0.6, 0.8, 0.9, 0.985. Shaded areas 
correspond to the variation of ipv in the interval — 8° < ipv < 8°; the lower and upper curves of the shaded area 
correspond to destructive and constructive interferences of ^^^^^^^(O) and ^^^^^^^(0), respectively. 

The measurement of the /o(980) signal in the 'jtt^tt'^ reaction gives the branching ratio BR{(j) 7/0 (980)) — 
(3.5 ± O.Stli) ■ 10"'' @, in the analysis of -/t:^tt^ and -fn+ir- channels it was found BR{(j} 7/o(980)) = (2.90 ± 
0.21± 1.5) -10-4 the averaged value is given in fT^]: Bi?(0 ^ 7/o(980)) = (3.4±0.4) • lO"''. Our calculations of the 
7r°7r° spectrum with the conventional Flatte formula (Section 4) support the statement of about the presence of 
large systematic errors related to the procedure of extracting the /o(980) signal. In our estimation of the permissible 
interval for the mixing angle cp, we have used the averaged value given by [|l6| , with the inclusion of systematic errors 
of the order of those found in M: BR{(P 7/o(980)) = (3.4 ± 0.4+J j) • 10"^ . 

The calculated values of BRQ) 7/o(980)) agree with experimental data for \ip\ > 25°; larger values of the mixing 
angle, \(p\ > 55°, correspond to a more compact structure of /o(980), namely, -R^(,(g8o) ^ 10 GeV~^, while small mixing 
angles \ip\ ^ 25° are related to a loosely bound structure, -^^^(ggo) > 12 GeV~^ (recall that for the pion R^. ~ 10 
GeV-2). 

The evaluation of the radius of /o(980) was performed in on the basis of GAMS data where the t- 
dependence was measured in the process iT~p — > /o(980) n: these data favour the comparatively compact structure of 
the qq component in /o(980), namely, -^^^(930) = 6 ± 6 GeV~^. 



III. CALCULATION OF THE PARTIAL WIDTH FOR <P{W20) 



7^V 



In the tt'^tt" spectrum, the resonance /o(980) is seen as a peak on the edge of the spectrum. So it is rather instructive 
to calculate the dependence of BR{(f){1020) — > 77r°7r") on A/tttt by using the conventional Flatte formula p^ . 

The resonance /o(980) in hadronic reactions is well described by the Flatte formula which takes into account two 
decay channels, /o(980) nn and /o(980) — > KK; the parameters of the Flatte formula may be found in P,p5[. 

In the calculation of the nn spectrum, we take into account the resonance contribution and the background, B(M.^^). 
The partial width dr0_>^^^ reads: 



dT 



0— >77r7r 



= r 



0^7/0 (980) 



ML) F,i_,(M2j 



Ml 



Ml 



M^Aml-Mj) F2,_,(M|) 



■ 2 dMl^ 
n 



(39) 



where Af/=0.98 GeV. Fthresh{M^^) is a factor which provides the threshold behaviour of the decay amplitude 



<^(1020) ^ 77r07r0 at M^ 



li. We have chosen it as FA 



is a scale parameter; below we put /iq — 2m.,^. 

Here we separate r^^^^^^^ggo) and introduce the ratio of phase spaces for [nn) 



threshiMD ^ 1 - exp[-(M2^ - ml)'^ / n^], where /zg 



7 and /o(980) + 7: 



t^^/o (980)7 



r(™W 



M2) 



(40) 



The values 5^ and are coupling constants for the transitions /o(980) 
spaces for the nn and KK states: 



nn, KK. Also p^rTr and are phase 



= ^/M^~^4^/{mnM^^) and p^R = V^^tt - ^m\/{l&nM^.^)- recall that 



in the case of M^^ < 4m|. the kaon phase space is imaginary, for a/M^ 



4m|- 



j-\/4rr? 



K 



M2 
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According to ||^,|2^, the Flatte formula is parametrized as follows: 



glp^^ ^ (0.12GeV)VM2,-4m2, g^p^^ ^ (0.27 GeV)^M2, - , (41) 

Mo = 0.99 ± 0.01 GeV. 

We parametrize the background term as B{M^^) = a + bM^^ where a and b are complex constants. The unitarity 
condition in the tttt channel tells us that the decay amplitude 0(1020) jirir has the complex phase factor related 
to the TTTT phase shift in the {IJ^'~' = 00++) channel: A^-^Tr = \A^-KTr\ exp [i(5o(M7r7r)]. Therefore, one has 



TTTT V 



+ S(M2^) ) e'^*o(Ai,,)_ (42) 



- Af 2^ - iglp^^ - igiiPKK 

We parametrize the phase A(5o(M7r,r) as /S.5q{Mt^t^) = Aq + Ai^AI-^T^/mo — 1) + A2(Af7r^/mo — 1)^ with nio= 1 GeV. 
The complexity of B{M^^) is determined by the difference of phases in the terms for the /o(980) production and 
primary background contribution. 

To calculate the 77r°7r° spectrum, one should multiply ( |39| ) by the factor related to the tt^tt" channel dF^^-yTrOTrO = 
idr^^^,^^. Figure 4 demonstrates Si?(0 ^ 77r07r°) calculated using Eq. (|9|), with BR{4) -> 7/o(980)) = 3.4 • lO"-*, 
5.4- lO^"' and Mq = 0.98 GeV, 1.00 GeV. The following parameters (in GeV units) are used for the background term 
B{Ml^) and the phase A(5g(M^^): 

(a) BR = 3.4 X 10"* , Mq = 0.98 : a = -1.24 + iO.74 , b = 1.91 - iO.95 , 

Ao = 99° , Ai = 11° , A2 = 440° ; 

(6) BR = 3.4 X lO"'^ , Mo = 1.00 : a = -1.96 + i0.073 , b = 2.37 + il.81 , 

Ao = 106° , Ai = -136° , A2 = 341° ; 

(c) BR = 5.4 X 10""^ , Mo = 0.98 : a = -0.84 + i0.91 , b = 0.73 - i2.63 , 

Ao = 98° , Ai = 203° , A2 = 735° ; 

[d) BR 5.4 X 10""^ , Mo = 1.00 : a = -1.13 + iO.69 , b = 1.39 - il.l3 , 

Ao = 94° , Ai = 9° , A2 = 462° . (43) 

It is seen that the data jsl agree reasonably with the Flatte parametrization for hadronic reactions given by |^,^ . A 
good description of the tPtt" spectrum at i?i?((/) —* 7/o(980)) = 5.4 • lO"'' supports the statement of about large 
systematic errors in the determination of the 7/0 (980) signal. 

IV. RADIATIVE DECAY /o(980) ^ 77 

The amplitude for the transition /o(980) — > 7(9^)7(5'^) has the same spin structure as — > 7/o(980): 

A^,, = e'Af^^^^iq^) gj^^ (44) 

where gl:^q'^ = and 5^+9^ 0. 
The invariant amplitude reads: 

where TVc = 3 is number of colours, z'^f^''^^^ — 5\/2/9 and -Z^Iq!!^^ = 2/9. The spin factor S'|*'j^^^(s, s', Q^) jg similar 
to that for </) — > 7/0, namely: 

SJX,{S. s', q') - -2m (^-s + s' + 4m' - --^f_^a(., q')^ . (46) 
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In the limit ^ 0, the term which makes different 5'f*'j^^^(s, s' , q^) and S'^*!^^^^ (s, s', g^) vanishes in the integral 
(^), and we have for yl/Q^^-y(0) an expression similar to (^8|): 



^/o-77(0) 



mJNrZ 



47r 



— V'/o(s)V'7(s) 
TT 



y/s(s — Am?) — 2to^ In 



(47) 



The photon wave function was fomid in the analysis of the transition form factors 7r° — > 7((5^)7, rj — s- j{Q'^)j, and 
?]' ^ l{Q'^)j it is shown in Fig. 5. With this wave function we calculate yl^^"_^^^(0) and Aj^^^^(O); these 



amplitudes plotted versus -Rjj,(g8o) are shown in Fig. 2b. 



The partial width F 



/o^77 



is equal to: 



™'|%-77(0)I^ 



(48) 



where 



^/o-»77(0) = COSipA 



nn ) 
/o-»77 



(0) 



sm^4^^. 



(49) 



Figure 6 demonstrates the comparison of the calculated Tf^, 



at different ^^^(ggQ) and (p, with the data ^ 



It 



is possible to describe the experimental data (F^„(98o)^77 = 0.28j;o;i3 |2^) at positive mixing angles as well as at 
negative ones: 77° < if < 93° and (—54°) < (f < (—38°). The use of the radius restriction (2) makes the allowed 
interval of mixing angles slightly narrower, see (3). 



V. RADIATIVE DECAYS (?!)(1020) -> 777, 77)', 77r°, 7ao(980) 

The decays (/)(1020) 777, 777', 777°, 709 (980) do not provide us with direct information on the quark content of 
/o(980); still, it looks necessary to perform calculations and comparison with data in order to check the reliability of 
the method. In addition, the decay 0(1020) — *■ 77r'' allows us to evaluate the admixture of the nn component in the 
(j) meson; as we saw in Section 2, this admixture affects significantly the value F^^^^^^ggQ). 

The amplitude for the transition 4> jP, where P — 77, 7;', tt*', has the following structure: 

Af_,^ = eA^^^p{q^)e^yapPaqi3 ■ (50) 

Radiative decay amplitudes </> — > ^rj,^ri' may be calculated similarly to the decay amplitude — > 7/o(980), see 
Eq. (p^), with necessary substitutions of the wave functions ipfo ^ V'r;, V'r;'! V'tt as well as charge and spin factors: 



^0-7P(O) = / — V'0(s)V'p(s)ln— , ^ „ ■ (51) 

47r J TT vs — V s — 4771^ 

The charge factors for the considered radiative decays are as follows: for the ss component in the reactions cj) — *■ 777, 777', 
■^^!^7^ = Z^^^^^n' = -2/3, and for the tt" production, Z^_^^o = 1. 

For the transitions (j) ^ jrj and — > 777' we take into account the dominant ss component only: — sin 0ss in 
77-meson and cos ss in 77'-meson, and sm9 = 0.6. 

The spin factors are equal to 

= S^^v'i^^ 9') = 47n„ s\ q') = 4m . (52) 

For pion wave function we have chosen — 2.0 GeV^^ that corresponds to i?^ = 10.1 GeV~^, the same radius is 
fixed for the nn component in 77 and 77'. As to the strange component in 77 and 7/, its slope is the same: b^j^ss) — 2 
GeV~^ that leads to a smaller radius, i?^(ss) = 8.3 GeV^^. 

The results of the calculations versus branching ratios given by PDG-compilation ||l^ are as follows: 

BR{^ 777) ^ 1.46 • 10^2 ^ BRpdg{<P ^ ?77) (1-30 ± 0.03) • 10"^ , 

Bi?(0 ^ 77'7) = 0.97 • 10-S BRpDG{4>^ri'^) = {{).&7tli\)-lQ-^ . (53) 
The calculated branching ratios agree reasonably with those given in |lq| . 
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For the process cj) jn'^ the compilation |^6| gives BR{(I) — > jir^) = (1.26 ± 0.10) • 10~^, and this value requires 
I sin(^y| ~ 0.07 (or \(pv\ — 4°), for with just this admixture of the nn component in (^(1020) we reach the agreement 
with data. However, in the estimation of the allowed regions for the mixing angle (p (see Fig. 3), we use 

|^y|=4°±4° (54) 

considering the accuracy inherent to the quark model to be comparable with the obtained small value of \(pv\- 

The process (/)(1020) — > 7ao(980) depends also on the mixing angle \fv\- the decay amplitude is proportional to 
sin</?y, namely, ^^^^^q^ = smipyA^^^^^. The amplitude A^"^^^^ is equal to that for the process 7/0, up to a 
numerical factor: A^^t.ao/^'^t^fo = Z^t^ao/Z^l^fo = ^ because Z;l^^^ = 1 (the value 3A^t^^^ is shown on Fig. 
2a). We have for the region RI^^q^^^ ~ 8 GeV-^ - 12 GeV-^: 

BR (0(1020) ^ 7ao(980)) ~ sin^ ipv ■ (14 ± 3) • 10"* (55) 

with lower values for ^^^^^j-ggo) ^ 8 GeV^^ and larger ones for ^^^^(ggo) ^ 12 GeV~^. At sin^ i^jy < 0.02, we predict 
Bi?(0(lO2O) ^ 700(980)) < 0.28 • lO''^. 

In recent paper |Q the partial width for 0(1020) jrj-K^ was measured, it was found i?i?((/)(1020) 'yrjiT^; M^i^ > 
900 MeV) = (0.46 ± 0.13) • 10"^. If at M^^ > 900 MeV the ratio of the resonance/background is of the order of the 
unity, that is quite possible, the value found in ||3l[| is in agreement with small value of \(pv\- 



VI. CONCLUSION 



Figure 7 demonstrates the ((/?, i?y^(ggQj)-plot where the allowed areas for the reactions 0(1020) — *■ 7/o980) and 
0(1020) 77 are shown. We see that the radiative decays 0(1020) — > 7/0 (980) and /o(980) — > 77 are well described 
in the framework of the hypothesis of the dominant qq structure of /o(980). For flavour wave function determined by 
^'/iaDour(/o(980)) = nn cos ((5 + ss sin ((5, the solutions found are either ip = —48° ±6° 01 ip = 85° ±5°. The solution with 
negative ip seems more preferable, that means the qq component is rather close to the flavour octet {ipoctet = —54.7°). 
Such a proximity to the octet may be related to the following scenario: the broad resonance /o(1530^25o)' which 
is, according to the if-matrix analysis, the descendant of the lightest scalar glueball after accumulation of widths 
of neighbouring resonances took the singlet component from a predecessor of /o(980), thus turning /o(980) into an 
almost pure octet. 

The dominance of qq component does not exclude the existence of the other components in /o (980) . The location of 
the resonance pole near the KK threshold definitely points on a certain admixture of the long-range KK component 
in /o(980). To investigate this admixture the precise measurements of the KK spectra in the interval 1000— 1150 
MeV are necessary: only these spectra could shed the light on the role of the long-range KK component in /o(980). 

The existance of the long-range KK component or that of gluonium in the /o(980) results in a decrease of the ss 
fraction in the qq component: for example, if the long-range KK (or gluonium) admixture is of the order of 15%, the 
data require either ip = -45° ± 6° or = 83° ± 4°. 

We thank D. V. Bugg, L. G. Dakhno, V. S. Fadin, L. Montanet and A. V. Sarantsev for useful discussions and 
illmiiinativc information. The work is supported by the RFBR grant 01-02-17861. 
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a b 

FIG. 1. a) Diagrammatic representation of the transition i^(1020) — > 7/o(980). b) Three-point quark diagram: dashed Unes 
I and II mark the two cuttings in the double spectral representation. 
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FIG. 2. Amplitudes for strange and non-strange components, ss and nn, as functions of the /o(980)-meson radius squared: 
a) <_:),jO)/4"4,„ and A^l%,^{^)IZ^%^, b) A'-;%{^) I Z^;% and ^-^^^(O)/^!:^^^. 
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FIG. 3. Branching ratio BR{(j){1020) 7/0(980)) as a function of the radius squared of the nn component in /o(980). 
The band with vertical shading stands for the experimental magnitude; five other bands, with skew shading, correspond to 
\>fi\ = 24°, 37°, 53°, 64°, 80° at -8° < ifiv < 8° . 
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FIG. 4. Branching ratios Bi?(</)(1020) — > 7r°7r°7) calculated with the Flatte formula (26) and parametrization given by 
(28) (solid curves) at different Bi?((/I)(1020) — » 7/o(980)) and Mf^(^^Q^. Dashed line shows the background contribution which 
corresponds to the elimination of the resonance term in (39) . The data are taken from [7] . 
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FIG. 5. Photon wave function for non-strange quarks, ip-y^nn{k^) = g-y{k^)/{k^ + rn^), where = s/A — m^; the wave 
function for the ss component is equal to ip-y^ssik'^) = g-i{k^)/{k'^ + fnt)'i the constituent quark masses are m=350 MeV and 
m»=500 MeV. 
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FIG. 6. Partial width V experimental data are from [22] (shaded area), a) Curves are calculated for positive 
mixing angles f = 77°, 80°, 85°, 90°, 93° and b) negative angles ip = -38°, -43°, -47°, -54°. 
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FIG. 7. The {(p, 7i^^(ggo))-plot: the shaded areas are the allowed ones for the reactions (?i(1020) — > 7/o(980) and /o(980) — > 77. 
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